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ABSTRACT 

Context. Planet formation and migration in accretion discs is a very active topic. Among the many aspects related to that question, dead 
zones are of particular importance as they can influence both the formation and the migration of planetary embryos. The ionisation 
level in the disc is the key element in determining the existence and the location of the dead zone. This has been studied either within 
the Standard Accretion Disc (SAD) framework or using parametric discs. 

Aims. In this paper, we extend this study to the case of Jet Emitting Discs (JED), the structure of which strongly differ from SADs 
because of the new energy balance and angular momentum extraction imposed by the jets. 

Methods. We make use of the (r, z) density distributions provided by self-similar accretion-ejection models, along with the JED 
thermal structure derived in a previous paper, to create maps of the ionisation structure of JEDs. We compare the ionisation rates we 
obtain to the critical value required to trigger the magneto-rotational instability. 

Results. It is found that JEDs have a much higher ionisation degree than SADs which renders very unlikely the presence of a dead 
zone in these discs. 

Conclusions. As JEDs are believed to occupy the inner regions of accretion discs, the extension of the dead zones published in the 
literature should be re-considered for systems in which a jet is present. Moreover, since JEDs require large scale magnetic fields close 
to equipartition, our findings raise again the question of magnetic field advection in circumstellar accretion discs. 

Key words, accretion, accretion discs — ISM: jets and outflows — Stars: formation. X-rays 



1. Introduction 

A detailed knowledge of the dynamics and chemistry of pro- 
tostellar accretion discs is unavoidable if one is to understand 
the initial conditions of planet formation and the environment in 
which their migration proceeds. The theory of protostellar ac- 
cr etion discs has been widely developed since the seminal work 
of lShakura & Sunvaevl ( ll973h and is commonly set in the frame- 
work of the standard accretion disc (SAD) in which viscosity is 
the sole agent to remove angular momentum and allow accre- 
tion (e.g. lPringld[T98l lD^ Alessio et al. l ll998l:lHueso&Guillotl 
l2005h . The current paradigm is that this is done through MHD 
turbulence, triggered and sustained i n the disc by the mag neto- 
rotational instability (hereafter MRI. lBalbus & Hawlevl[T99 1). 
This instability requires a good coupling between the disc 
, material (mostly neutrals) and the magnetic field, therefore re- 
quires a minimum degree of ionisation to exist. If the ionisa- 
tion level is not high enough for the MRI to exist, the viscos- 
ity is very low: this results in a very low level of accretion in 
the correspondi ng region of the disc, then termed "dead zone" 
jGammidl 19961) . Dead zones may play an important role in the 
field of planet formation and migration. Rossby vortices can be 
triggered at the interface between an active and a dead zone, 
trapping solids within them and leading to a burst o f planet for- 
mation dVarniere & Taggeil 120061: iLvraet al.ll2009l) . They also 
have been put forward as an agent able to halt planetary mi- 
gration, a long-standing issue in explaining planetary systems 
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such as our own ( iMasset et al.l 120061: iMatsumura et alj|2009l) . 
Locating dead zones in an accretion disc is therefore an im- 
portant task. Several authors have studied this question, gener- 
ally using either i) a passive disc (Glassgoldet al. ( 1997) uses a 
powe r-law disc, [M atsumura & Pudrit j (120031) the passive disc 
from Chiang & Goldreicb (J 9971) ). or ii) a standard accretion 
disc ( iFromang et alj 20021) . These studies show that a dead zone 
can generally exist in the inner parts of accretion discs, typically 
extending from 0. 1-1 AU up to a few AU. 



Another way to indirectly probe the physics of the inner disc 
regions is to look at the jets from these sources. Bipolar jets are 
fundamental to star formation and are observed from the ear- 
liest (Class 0) to the latest stages (classical TTauri stars, Class 
2) of the process. They are believed to originate from the in- 
ner disc (up to a few AU), precisely where the aforementioned 
works predict the existence of a dead zone. Extended MHD disc 
winds are, to date, t he best candidates to explain jet kinematics 
dFerreira et alj2006h . Because they change the energy balance of 
the disc, the jet emitting disc (JED) structure differs from that of 
the SAD one. This was shown in Combet & Ferreira ( 2008) (pa- 
per I hereafter) and questions the existence dead zones in the in- 
ner parts of jet launching discs. This is the issue we address in the 
present paper which is outlined as follows. We start by a quick 
recap of disc wind solutions and JED structures in section |2] 
Section[3]presents the procedure we follow to calculate the ion- 
isation level of a JED while the results are shown in section |4] 
Finally, discussion and conclusions are drawn in section|5] 
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2. Jet Emitting Discs: a quick reminder 

It is widely accepted that protostellar jets observed in young 
stellar objects (YSOs) are launched from the surface of accre- 
tion discs via the mag neto-centrifugal launching mechanism, 
originally developed bv lBlandford & Pavnel(ll982h . Solutions to 
the problem have gen erally been found by treating the disc as 
a bou ndary condition jBlandford & Pavndl 19821 : IVlahakis et al.l 
I2OOOI) . hence preventing any quantification of the interplay 
between the disc and the jet. This difficulty has been lifted 
in the self-similar framework of the Magnetised Accretion- 
Ejection Structures (MAES) model, where the set of steady- 
state MHD equations are consistently solved, going from the 
resis tive regime withi n th e disc to t he ide a l MHD wind/jet 
jFerreira & P elletier 1991 iFerreiral 119971: ICasse & Ferreiral 
l2000al:lFerreira & Cass3l2004 . MAES rely on a self-sustained 



turbulence to provide the required anomalous transport of mag- 
netic field (through an effective magnetic diffusivity), along with 
magnetic fields smaller than, but close to, equipartition (i.e., the 
plasma /3 = 2iJ.oP/B~ > 1). This may appear problematic as 
the magneto-rotational instability, invoked as the source of tur- 
bulence, is known to be quenched whenever the field reaches 
equipartition. This difficulty is lifted when remembering that (1) 
MAES are cl ose to but always below t he marginal stabihty hmit 
for the MRI dBalbus & Hawlevll 19911: iGammie & Balbuslll994 
(2) th e MRI could still persist at equipartition ( Blokland et al.l 
I2OOI . Note also that other sources of MHD turbulence could 
be at work around equi partition in a medium as complex as a 
stratified accretion disc (iKeppens et al.ll2002h . 
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Fig. 1. Log-density (cm^^, colour contours) and streamlines 
(solid black lines) of a JED powering jets that seem consistent 
with observations (see text for more details). The dashed yellow 
line shows the position of the slow magnetosonic (SM) point (lo- 
cated at zsM = 1.67 h). The SM point occurs in the ideal MHD 
regime above the disc, whereas the flow is deflected vertically 
(ie Ur - 0) below, in the underlying resistive layers. Note that in 
order to resolve the disc internal structure, the vertical scale has 
been expanded. 



2.1. Self-similar MAES solutions 

In the present study we use a disc wind solution unde rgoing a 
surface heating described in ICasse & Ferreiral (l2000bl) . Due to 
this heating, the vertical plasma pressure is increased which en- 
hances substantially the jet mass load. This additional heating 
was assumed to be powered by dissipation of MHD disc tur- 
bulence and parametrised as a tiny fraction of the disc accre- 
tion power (/ = 8 X 10 in the solution chose n here). This 
new c lass of solutions, referred to as "warm" in iPesenti et al.l 
(12004 . is still "cold" in a dynamical sense as the initial thermal 
energy remains negligible compared to gravity, and wind accel- 
eration is still mostly magnetic. The disc ejection efficiency ^ is 
defined by Mace f^. While isothermal disc wind models pro- 
vide ^ around 0.01, these warm models can achieve larger val- 
ues. For the solution used here, we have ^ - 0.04 from a disc 
of aspect ratio e - h/r = 0.03, where h{r) is the disc half thick- 
ness. If the JED is established between say, ri„ - 0.04 AU and 
rj - 0.3 - 1 ALQ this translates into an ejection to accretion 
mass rates ratio 2Mjet/Macc - ^In ^ = 0.08 to 0.13, in agree- 
ment with observations. The asymptotic poloidal jet velocity is 
a factor V^l-S larger than the keplerian velocity at the field 
line footpoint, where /I ^ 1 + 1 /2^ is a measure of the magnetic 
lever arm ([perr eira 1992). With A - 13.15, this model is found 
to reproduce well several observed properties of atomic T Tauri 
jets such as the rotation signatures, the ejection to accretion ra- 
tio, and t he centroid velocities in atomic spectra along the jet 
axis (e.g. IPesenti et al.ll2004l:ICabritll2007l) . 

' The question of the radial extension of the launching region of 
the disc has been debated at length. If X-wind confi gurations, wher e 
the jet comes from the innermost region of the disc iShu et aLll2000l) . 
could surely contribute to the mass loss, an extended disc-wind com- 

i ionent appears indis pensable, especially in the younger Class I objects 
Perreira et alj2006h . 



The resulting (r,z) distribution of density is shown in fig.[T] 
along with some streamlines (black solid lines). Note that the 
range of launch radii produces an "onion-like" structure with a 
fast, dense axial beam surrounded by progressively slower and 
wider streamlines that recollimate on larger scale, producing a 
high apparent jet collimation. This is however not shown here as 
this work focuses on the properties of the Jet Emitting Disc only. 
As can be seen on the figure, only the upper layers of the disc 
are deviated into the outflow, while most of the material ends up 
at the inner boundary of the disc, to be accreted by the protostar 

The global picture of a protostellar accretion disc would then 
be the following (see fig.l of paper I): i) in the inner parts from 
(rin to rj), a jet emitting disc (JED) having MAES properties and 
ii) a standard accretion disc (SAD) from rj to rout- As mentioned 
above, SADs rely solely on the radial viscous transport of an- 
gular momentum for accretion to proceed. Conversely, MAES 
solutions show that as soon as a powerful jet is present, most 
of the angular momentum and energy within the disc are evac- 
uated vertically in jets. This suggests that JEDs should have a 
very different structure than SADs and is developed hereafter. 

2.2. The jet emitting disc's structure 

In paper I, we analytically derived the radial properties of JEDs 
for a given protostellar mass M and mass accretion rate in the 
disc Mace- This was done in a formalism analogous to that of 
SADs, but using the energetic balance suggested by the MAES 
solutions. It was found that for a given M - Mace combination, 
a JED is cooler, lighter and thinner than its equivalent SAD (we 
refer the reader to that paper for more details). 

Its lower temperature should allow for the presence of a JED 
to be tested by the spectral energy distribution (SED) of the disc. 
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JED+SAD, z^^=OJ R ., 
SAD, z^^ = 0.7 R 
SAD, i^^ = 
JED+SAD, z = 




Fig. 2. SED for different disc configurations with M^cc - 
10"^ Mq yr Solid lines correspond to a disc composed of a 
JED up to ry = 0.3 AU and a SAD in the outer regions. Dashed 
line represent the case where the SAD fills all the disc, up to the 
innermost radius = 0.04 AU. For each case, the source of ir- 
radiation (UV from accretion column hot spot) had been located 
in z = (thick lines) and z = 0.7/?star (thin lines). 



This was illustrated using very crude assumption^ for the disc 
geometry in paper I. Therefore, we present an updated version of 
the SED in fig.|2]for illustrative purposes. Both viscous heating 
and illumination from the central object are included in the SED 
calculation. Assuming most of the irradiation comes from the hot 
spot of the accretion column onto the protostar, we calculate the 
SED for two altitudes of the hot spot, Zsrc = and Zstc - 0.7/?*. 
We consider both the cases where the SAD extends down to ri„ 
and where it is supplemented by a JED in the inner parts, with 
rj - 0.3 AU. As can be seen in fig. [31 for a given spot altitude, 
there is a clear signature of the presence of the JED. However, 
it appears very hard to distinguish between a JED+SAD con- 
figuration with Zsrc = O.lRi, and a single SAD with Zsrc - 0. 
This simple example illustrates that one cannot solely rely on 
the SED to test for the presence of a JED in a given system. 
But this indirect method should nevertheless prove useful when 
used in combinatio n with other diagnostics (e.g. jet kinematics, 
iFerreira et al ] l2006h . 

3. Calculating the disc ionisation 

Several agents can contribute to the ionisation of an accretion 
disc. In his seminal work on dead zones. iGammiei (i 1 996.) consid- 
ered ionisation by the cosmic ray background only, of a pow er- 
law disc. Using a passive disc. lMatsumura & Pudritj (l2003h in- 
cluded in their calculations ionisation from cosmic rays. X-rays 
and radioactive decays (that remains marginal) in the disc. They 
concluded that cosmic ray (CRs) ionisation was the leading pro- 
cess while the X-ray contribution could become dominant if the 
X-ray photon energy was sufficient (5 - 10 keV). This result 
certainly holds if CRs can reach the disc. However it has been 



^ In fig. 4 of paper I we calculated the SED assuming a given incident 
angle for the radiation. This lead to an overestimation of the irradiation 
flux in the inner parts of the disc making it too bright at small wave- 
length. This is corrected here by using the calculated JED geometry 
consistently to estimate the incident angle at each radius. The qualita- 
tive conclusions of paper I are unchanged, simply the actual values of 
the flux. 



argued that they may not be able to do so because of the winds 
associated to star formation, in which case. X-ray ionisation be- 
comes fundamental again (Glassgold et al. 1997; Fromang et al] 
This is particularly true for JEDs that are, by definition, 
the source of a strong wind that eventually recollimates into a jet. 
Therefore, we only focus in this paper on the X-ray ionisation of 
the JEDs. 



3.1. Ionisation rate from X-rays 

A well known result from X-ray ionisation studies is that the 
ionisation rate is dominated by secondary ionisation, i.e. the ion 
pairs co llisionnaly produ ced by a fast primary photo-electron 
( Krolik & KallmanI 1 9831). The X-ray ionisation rate is written 
as (iMatsumura & Pudritzll2003l and references therein) 



fx 



AncfkTx 
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cr{kTx) 



X " exp 
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(2) 



where Lx is the source X-ray luminosity, d is the distance from 
the source, kTx the energy of the photons, <T(kTx) is the absorp- 
tion cross section, Ae ^ 37 eV is the average energy required for 
ion pair productior(3 and J(xo, t) is the attenuation factor The 
latter is calculated using x = E/kTx and is function of the op- 
tical depth, T(kTx) = NncrikTx) with Nh the column density 
along the line of sight, and of the cut-off energy xq = Eo/kTx. 
The absorption cross section is energy dependent and expressed 



o-{kTx) = o-Q 



kTx V 
Ikevj 



with o-Q ^ 8.5 X 10-^^ cm^ and n ^ 2.81 dGlassgold et alj|1997l: 
llgea & Glassgoldll999l:lF^omang et al.ir2002h . 



3.2. Recombination processes 

Electrons in an accretion disk may recombinate via several chan- 
nels, by interacting with both the gas and dust grain phases. 
Dissociative recombinations in the gas phase have been exten- 
sively studied both theoretically and experimentally, giving a 
relatively well-defined picture of recombination with molecular 
or me tallic ions (for a review, see iFlorescu-Mitchell & Mitchell 
I2006h . 

The picture becomes however far more complicated when 
gas-grain chemistry is considered. Several authors have investi- 
gated the importance that recombinations of electrons on the sur- 
face o f dust grain s may play on t he ionisation l evel of accretion 
discs (ISano et al.l|2 000; Semen ov et alj 120041: lllgner & NelsonI 
120061: IWardlell2007t iSalmeron & Wardlel l2008>). It is generallv 
found that including grain chemistry decreases the ionisation 
level of the disc (with respect to the results of gas-phase cal- 
culations only), provided that the grains are small enough (< 
1 urn ) and well- mixed with the gas, on the whole disc thick- 
ness ISano et alj 120001: lllgner & NelsonI |2006|) . These condi- 
tions may not always be fulfilled as theoretical considerations 



^ The parenthesis should read (kTx - iP)/^^, as the energy of the 
primary photo-electron depends on the ionisation potential (IP) of the 
medium. However, /P is ~ a few eV whereas kT^ > 1 keV. 
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and SEDs of TTauri accretion disks suggest that grains grow 
and settle onto the disc midplane during the protostellar evolu- 
tion (Dullemond & Dominik 2004; D'Alessio et al. 2 006) Note 
also that in most of the models in illgner & Nelson! (l2006h . the 
gas-only and gas-grain chemistry converge in the inner disc 
(< 0.6 AU), where dust grains may not survived. Finally, the 
detailed numerical studies investigating the effect of dust grains 
on the MRI tSalmeron & Wa rdle 2008; Turner et al. 2010), gen- 
erally do so at distances of a few AU that are not relevant to the 
JED. 

Given the complexity and uncertainty surrounding the grain 
chemistry and distribution in the discs, we restrict ourselves to 
recombination processes in the gas-phase for the remainder of 
this paper. This choice is motivated by the fact that JED are rele- 
vant only to the innermost regions of the discs, where dust grains 
may not survive and where it is difficult to extrapolate the results 
of existing studies. 

3.3. lonisation fraction 

The ionisation fraction - tie/nn is generically obtained from 
the solution of the evolution equation of the electronic density «g. 
Considering the gaseous phase only, this evolution is des cribed 
by system of coupled equations (e.g. lFromang et aill2002h 



drie 
at 



(3) 



IT 



where hm and «m+ are the density of metals and metal ions re- 
spectively, «„,+ is the density of molecular ions, the ionisation 
rate is given by eq. ([T]i, /3 is the dissociative recombination 
rate, /3r is the radiative recombination rate and yS, is the rate of 
charge transfer between molecular ions and metal atoms. 

Steady-state, that is assumed in a majority of the studies pub- 
lished in the subject, is valid if the ionisation and recombination 
timescales are shorter than the accretion timescale in the disc. 
This is true for SADs as they are characterised by a slow inward 
radial velocity, hence long accretion timescales. Also assuming 
no metals («m = «m+ = 0, rig = «m+), eq- ® leads to 




(4) 



A priori, more caution is necessary in the case of JEDs where 
the accretion proceeds much faster, namely at sonic velocities. 
For a JED, the accretion timescale is therefore given by 



— {Cs is the isothermal sound speed) 



(5) 



while the ionisation and recombination timescales obtained from 
eq. (O read 



^lon — ^ 

and 

1 



PXeTlH 



(6) 



(7) 



In steady-state, these two latter timescales are equal, fion = free- 
For our calculations, we first assume that steady-state holds and 
we check a posteriori (in ^4.11 ) the extent to which this hypoth- 
esis is valid. 



3.4. Dead zones and critical ionisation rate 

The MRI, that arg uably provides the disc's a nomalous viscos- 
ity and resistivity (iLesur & Longarettill2009t iGuan & Gammid 
2009), only sets in if the Alfvenic time scale ta = h/VA is shorter 
than the diffusion time scale - h^/rj due to electron-ions col- 
lisions. Given the Alfven velocity of the medium Va and the disc 
scale height h, this criterion translates into a magnetic Reynolds 
number = hVA/rj that must be greater than a critical value 
The value of the critical Reynolds number is known to de- 
pend on various factors, such as the presence of a net magnetic 
flux through the disc, but is not firmly established yet. Therefore, 
we adopt the approach widely found in the literature and con- 
sider both 'K™' - 1 and '??™' = 100. Using the Ohmic resistivity 
(Blaes & Balb'uslll994h 

T] = 234 cm^s"' , 

Xe 

and th e fact that a JED requires a magnetic field close to equipar- 
tition dFerreira & PeUetieii[T99l . this condition becomes! 



(oil) (t^Is*?; 



(8) 



where e - h/r is the disc aspect ratio. From eq. (|4]i and us- 
ing the disc vertical hydrostatic equilibrium, the definition of 
the accretion rate and the c oUisional recombinat ion rate y6 = 
3 X 10 'T '/^ cm^ s ' (e.g. iFromang et al.ll2002h one gets the 
ionisation fraction to be 



M 



1/2 



M 



-1/2 



llAU/ 



,10-^Moyr-i, 
Inverting eq. dS), this leads to a critical ionisation rate 



(9) 



M 



lO-^Moyr-i , 



(lAu) 



(10) 



Note that this result is obtained for the typical recombination rate 
/3 mentioned above. 

3.5. X-ray properties of young steiiar objects 

Young stellar objects of all stages (Class to III) are known 
to be X-ray emitters. In the 1980's, the first discoveries were 
made by the Einstein X-ray ob servatory in the less embed - 
ded objects (Class II and III, e.g. lFeigelson & Decamplilll98lh . 
The ASCA and ROSAT sateUites aUowed the detection of 
X-ray s from Class I pro t ostars during the f ollowing decade 
(e.g.. iKovama et al.l [19961 iGrosso et al]ll997h . More recently, 
the Chandra Orion ultra-deep project (COUP) survey, re- 
vealed X-ray detectio n in the youngest Class objects (e.g. 
iPrisinzano et al.ll2008h . 

YSOs typical X-ray luminosities range between 10^** - 
10^" erg s"', with flares reaching 10^^ erg s"' and with typical 
photon temperatures of kTx = 1-5 keV. Quiescence properties 



jammie 



ivjaiiiiiiig II177UII gives an equivalent expression of %„ for a SAD 
using iHavashi iigsih 's resistivity and the central disc temperature T„ 



instead of the disc aspect ratio e. 
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may vary slightly from one class of object to another but typi- 
cal values are Ly™^'' = 10^^ erg s~' and ^r^""'^'^ - 1- 3 keV 
(llmanishi et al.ll2003l: iFavata et al.ll2005t IWolk et al.ll2005l). 

Using the statistics of the COUP survev. l^lk et al.1 ( l2005h 
determined a rate of one flare per week per star, with flare du- 
rations ranging from a few hours up to 3 days. The flare lu- 
minosities, spectral properties and variability distributions ap- 
pear fairly constant through all the stages of protostella r evolu- 
tion, suggesting single underlying mechanism (Prisinza no et al.l 
|2008|) . Fiducial flare characteristics are L^""*^ - 10^^"""" erg s"' 
and a spectrum harder than in quiescence, kT^^'^^ > 5 keV. 

One might expect these flares to be due to reconnection 
events in magnetic loops connecting the proto star to its cir- 
cumstellar disc. However, iGetman et al.l (l2008bh found no dif- 
ferences between accreting and disc-free systems (i.e. Class II 
and Class III). Thus, these reconnection events must be related 
to magnetic loops in the stellar magnetosphere alone. The size 
of these loops, hence the location of the X-ray emitter, can be 
estimated from the flare propertie s but remain "flare-mo del" de- 
pendent. From the COUP survey, iGetman et al.l (l2008ah obtains 
half-lengths of L ~ 10'^ - 10'^ cm. 

From the above, we see that the X-ray source alone intro- 
duces three parameters in the problem: i) its luminosity Lx, 
ii) the photon temperature kTx and iii) its location/altitude Zsrc, 
linked to the magnetic tube length. These three parameters will 
be systematically varied within their (reasonable) respective 
ranges. 

3.6. Choosing the appropriate self-similar solution 

Using eq. ([T]), we wish to calculate the ionisation rate of a JED 
for a given set of M - Mace- To that end, the appropriate MAES 
solution needs to be selected. On the one hand, each M - M^cc 
gives rise to a unique JED aspect ratio profile e(r), as shown in 
paper I (see also fig. |8]l. Such a profile is in particular strongly 
dependent on the dominant opacity regime, which varies along 
the radius. On the other hand, a self-similar MAES solution is 
characterised by a specific value of e. Therefore, for a given sys- 
tem (M - Mace), we have to select a self-similar solution whose 
e best matches that of the JED radial structur^ 

Our system's canonical physical parameters are inspired 
from the Class I/II object DG Tau and read M - 0.5 Mq and 
Mace - 10"^ Mq yr For these values the JED radial profile 
calculation gives an average value over the JED extension of 
e ~ 0.03 (Paper I). As described in Paper I, this calculation is 
valid in the optically thick regime which extends, for these par- 
ticular parameters, up to r ~ 1.8 AU. The MAES solution shown 
in fig.[T]is characterised by e = 0.03 and is the reference solution 
that will be considered in the following section. We also arbitrar- 
ily choose a radial extension of ry = 1 AU only, well inside the 
validity region of the JED calculation. 



4. Results 

Before exploring the question of the existence or absence of a 
dead zone in JEDs 04. 2b . we first present the validity range of 
the steady-state hypothesis we make. 



^ There naturally is a degeneracy in M - Mace, several combinations 
of which giving rise to approximately the same values of e. However, 
this is not a point we wish to explore in the present paper. 




R(AU) 



Fig. 3. Accretion (dashed) and ionisation (solid) timescales on 
the disc midplane as function of distance from the central ob- 
ject, with the X-ray source {Lx = 10^^ erg s"') located on the 
midplane (zsrc - 0). Two electron temperatures have been con- 
sidered. Calculations done for our reference MAES solution. 



4. 1 . Timescales and steady-state 

To infer the ionisation fraction from the ionisation rate using 
eq. (21, steady-state is assumed. Once the ionisation rate and 
fraction are known, we can verify a posteriori whether or not this 
hypothesis holds (fion < face), using eq. (|5]l and eq. (|6|l. Although 
this check can only be done once the calculations in ^4.2l are per- 
formed, it seems more relevant to present these considerations 
beforehand. 

As evoked by fig. [1] the problem should be treated in a full 
2D configuration. However, we present several results at the disc 
midplane. Indeed, it is the hardest part of the disc to ionise be- 
cause of the amount of material the X-rays have to cross, and is 
therefore the most constraining location to look at. 

The radial evolution of the accretion and midplane ionisa- 
tion timescales are plotted in fig. [3] When the system's physical 
properties are fixed (M, M, e), the accretion time (dashed line) 
simply scales as face r^^^- The solid lines represent the evolu- 
tion of the ionisation timescale, for a source located in Zsic = 
with Lx - 10^^ erg s '. For a 3 keV photon (thick solid line), 
the ionisation time is always shorter than the dynamical one 
which validates the use of the steady-state for this configura- 
tion. This is not the case for the 1 .5 keV photon (thin solid line) 
where the ionisation time is longer than the accretion time when 
r > 0.1 AU. This is just an example illustrating a situation where 
the full time-dependent treatment should be used. 

Figure |4] summarises the validity regime of the steady-state 
hypothesis at r = 1 AU (again, more constraining than any 
smaller radii). In the luminosity-photon energy plane, the solid 
black and blue lines locate where face = hon for Zsrc - and 
Zsrc = 3/Jo respectively. Above the line, steady-state can be as- 
sumed whereas it fails below it. The frontier between the two 
regimes naturally changes with both luminosity and photon en- 
ergy: a high luminosity compensate for a small photon temper- 
ature and vice versa. The source altitude is also a critical pa- 
rameter: most of the parameter space is covered by the steady- 
state hypothesis for Zsrc = 3/?o (blue line), whereas it is not 
the case for Zsic = (black line). For the latter, a time depen- 
dent approach should be used for the lowest photon temperatures 
{kTx < 2 keV). 
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Fig. 4. Black and blue solid lines represent face = fion in the 
luminosity-temperature plane for an X-ray source located at 
Zsrc - and Zsrc - ^ Rq respectively. Steady-state is valid only 
if 4cc > fion, namely above each line. The shaded area provides 
the most conservative estimate. Calculations are done with our 
reference MAES solution. 



The validity of the steady-state hypothesis depends on the 
X-ray source parameters (luminosity, energy, location). Overall, 
this approach holds on most of the parameter space. In the fol- 
lowing we restrict ourselves to steady-state calculations, keeping 
in mind they break down for unfavourable combinations of the 
parameters. 

4.2. lonisation ofJEDs 
4.2.1. lonisation rates 

The ionisation rate of the MAES solution presented in fig. [T]is 
calculated using eq. ([T). A 500x500 MAES box ensures the con- 
vergence of the column density numerical integration necessary 
for the optical depth calculation. The attenuation factor J(xo, t) 
given by eq. Q is computed from xp = IkeV /kTx to.Xmax - 100 
jMatsumura & Pudrit j2003h . 

Maps of the log-ionisation rate of our reference MAES struc- 
ture are given in fig.|5]for a X-ray source located in Zsic - (left 
panel) and z^ic = ^©(right panel). As for the density map, the 
black area near the vertical axis is the region where the self- 
similar solution is not defined. The ionisation rate spans over ten 
orders of magnitude, from ~ 10""^ s"' in the wind to ~ 10"^" s"' 
in the disc. The importance of the X-ray source location is vis- 
ible when comparing the two maps. Moving slightly the source 
from Zsrc = to Zsrc - Rq is enough to increase the ionisation 
rate of the disc by ~ two orders of magnitude. As will be seen 
later, going to higher altitudes (reconnection of magnetic loops) 
has an even more drastic effect. 



4.2.2. Dead zones 

In this section, we compare the disc midplane ionisation rate to 
the critical rate given by eq. ( fTOb . Note that our aim in doing 
so is not to determine with precision the geometry of the dead 
zone (thickness, extension) but rather to hint at the possibility 
for a JED to host one. For that reason, all the results we present 
here are located at the disc midplane, which is the hardest region 
of the disc to ionise. If (xlCx^ ^ 1' mid-plane is MRI-active 



and no dead zone is present anywhere, whereas the opposite case 
tells us that the disc may harbour a dead zone, without specifying 
its thickness. 

The canonical parameter configuration is Lx - 10^^ erg s"', 
kTx - 3 keV and - 1 . We systematically explore the param- 
eter space and results are presented in fig. |6] Each panel shows 
the radial variation of fx /f™' ^^'^ horizontal line repre- 

sents the critical value f^/f™' ~ ^ (corresponding to - 1). 
Note that the critical rate does not depend on the X-ray source 
properties and remains constant (at a given radius) throughout: 
a change in the fx/f™' ratio translates directly into the same 
change in (x- 



X-ray source location — fig. [6] upper left panel As already 
briefly shown in fig. |5] the location of the source w.r.t. the disc 
midplane is an important parameter This is made clearer by the 
upper left panel of fig. |6] where the ionisation rate varies over > 
4 orders of magnitude when the source moves from Zsrc = to 
Zsic = 10 (reconnection in magnetic loop). A change in the 
radial profile is also noticeable due to the change of geometry 
with the source acting as a lamp post on the entire disc for the 
highest altitudes. For all these configurations, the steady-state 
hypothesis is valid (see fig.© and the disc appears slightly sub- 
critical (i.e. with a possibility of dead zone) only when Zsm - 0. 
But this is inconsistent with the way we believe these X-rays are 
actually produced. 

Luminosity — fig.|6]upper right panel The ionisation rate is di- 
rectly proportional to the luminosity of the source and the effect 
is therefore straightforward. The source is located in Zsic = 
and for the lowest luminosity, a significant part of the disc 
should harbour a dead zone. However, for this low luminos- 
ity, the steady-state hypothesis is only marginally valid and a 
time-dependent approach should be used to confirmed this re- 
sult. Nevertheless, even for the lowest luminosity used, the dead 
zone disappears as soon as the source is displaced at a slightly 
higher altitude (not shown here). 

Photon temperature — fig.[6]lower panels The ionisation rate 
is also very sensitive to the photon temperature. The two bottom 
panels in fig.|6]show the evolution of the fx/f™' ^^^^'^ when kTx 
is changed between 1 to 5 keV. First of all, with a source in 
7,j.^ = 3Rq (left), the MRI is active in the entire JED if kTx > 
2 keV. For lower energies, the disc midplane appears strongly 
sub-critical, implying the existence of a dead zone. However, 
for this lowest temperature, the a posteriori verification showed 
that this was a regime where steady-state could not be assumed 
(fion > face) and wc cannot conclude in the presence of a dead 
zone without a full time-dependent study. 

This timing issue is solved when the X-ray source comes 
from some reconnection event (zsrc = 10 Rq, right panel). In that 
case, the ionisation rate is always larger than its critical value, 
whatever the energy of the X-ray photons. No dead zone can be 
present. 



Recombination process In order to calculate the ionisa- 
tion fraction Xg, we have assumed that dissociative re- 
combination with molecular ions was the main mechanism, 
and parametrised it with an effectiv e value ~ 3 x 
IQ-6J--1/2 ^,^3 g-i ^ Glassgoll 119991: iFromang et all 

120021 iMatsumura & Pudritz<i2003i) . This value actually depends 
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Fig. 5. Log-ionisation rate (s"') of the self-similar accretion-ejection solution presented in fig. [T] The system's characteristics are 
M - 0.56 Mq, M - Mq yr"', Lx - 10^^ erg s"' and kT., - 3 keV. Left: the X-ray source is located in z^rc - 0. Right: the X-ray 
source is located in Zsic - Rq- Black lines are iso-contours. 
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of the composition of the gas, each species having its own re- tion (e.g.. lSano et al ] l2000l:LSemenov et alJl2004l) . but this is be- 

combination rate. Theoretical and experimental studies give a yond the scope of the present paper When metals dominate, ra- 

range of individu al parametrisations (i - /jpT " (wi th n of- diative rather than coUisional recombina tion occurs, in which 

ten found ~ 0.5, iFlorescu-Mitchell & Mitchelll l2006l) . Using case P = Pr = 3 x l0-"r -'/2 cm^^ s ' (iFromang et al J 120021: 

these values, several authors have treated the full, or a re- iMatsumura & Pudritzll2003l) . 
duced, chemical reaction network in their dead zone calcula- 
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In order to evaluate the sensitivity of our results to the re- 
combination process, a simple exercise consists in using a para- 
metric expression for the recombination coefficient yS = /3oT'^^^ 
and changing the normalisation factor f3o from 3 x 10"" to 
3 X 10"^ cm"* s"' . This does not provide us with quantifica- 
tion for a realistic process but rather gives a feeling on how the 
situation changes with this coefficient- For this reason, we do 
not present a figure but simply summarise our findings below. 
For the canonical parameters Lx - 10^^ erg s"', kTx - 3 keV, it 
is found that: 

- For a source in Zsrc - 0, the steady-state hypothesis is valid 
on the entire radius range for /3o >3x 10^^ cur' s"' , but fails 
below. In this regime, no dead zone is found, whatever /3o. 

- If Zsic = 10 Re, there is no restriction to the steady-state 
assumption and the disc midplane is always MRI active. 

Of course, these conclusions are dependent on the X-ray lu- 
minosity, the photon energy and source location. Nevertheless, 
even for the lowest energies and luminosities, we do not find any 
dead zone present in the disc, whatever /3o, if the X-ray source 
comes from a reconnection event in Zsic = 10 Rq. 

Critical magnetic Reynolds number In ^3.41 it was noted that 
the value of the critical magnetic Reynolds number is still largely 
unknown. Figure|6]shows the results obtained for'??™' = 1. The 
critical ionisation rate given by eq. ( fTOl i scales with C??™')-. All 
other things being equal, fig.|6]can also be used to visualise the 
results obtained for '??™' - 100 simply by shifting the red "crit- 
ical" line from 1 to 100^ - 10'^. This configuration is clearly 
less favourable to the MRI. With Lx - 10^^ erg s only hard 
flares with kTx ^ 3 keV and located in Zsi-c 7/?© will pro- 
vide sufficient ionisation for the MRI to be present in the entire 
disc. However, these limits still accommodate most of the X-ray 
flare properties observed in young stellar objects. We therefore 
conclude that, even if ^^Jj"' = 100 makes it easier for a JED to 
host a dead zone, this nevertheless remains unprobable given the 
observational X-ray properties of YSOs. 

4.3. Dependency on the underlying MAES solution 

All previous results on the ionisation structure of JEDs have been 
obtained with our reference solution. Are they strongly depen- 
dent on the solution used? 

For the set of "w arm" (with additional heat deposition at 
the disc surface lavers lCasse & Fe rreira 2000a) and "cold" (e.g. 
isothermal magnetic surfaces Iperreirai 1 997.) solutions we have 
available, we calculate the disc midplane ionisation in the way 
described above. For a fixed stellar mass, accretion rate and X- 
ray source we find that: 

- For a given e, all warm solutions give the same disc midplane 
ionisation degree whatever ^. 

- The same applies to all the cold solutions. 

- For given e and ^, warm and cold solutions give similar re- 
sults. 

- For a given ^ but different e, the disc ionisation changes 
strongly. 

These results show that the wind properties, be it through 
the extra-heating term of the warm solutions or the value of the 

* Note that when the recombination coefficient changes, so does the 
critical ionisation rate which is then not given by eq. ilOl anymore. 
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Fig. 7. Solid black lines: contours of the optical depth t for a 
1 keV X-ray source located in Zsic = 0.01 AU. The values of r, 
from top to bottom are t = 1, 10, 100, 1000, 3000 and 5000. 
The colour scale represents the log-density (cm"^) of our refer- 
ence solution, as in fig.[r| 



ejection index ^, play a very little role on the ionisation state 
of the disc. We do not show the corresponding plots as the al- 
most superimposing lines give no more insight that what has al- 
ready been said. The equivalence between cold and warm so- 
lutions w.r.t to the ionisation of the disc is easily understood 
as the wind is very tenuous, no matter what the configuration, 
and brings therefore a only a small contribution to the column 
density between the source and the disc. We illustrate the lat- 
ter point in fig. Q where the solid contours represent the op- 
tical depth T{kTx = 1 keV), for an X-ray source located in 
Zsic - 3/?0 ~ 0.01 AU, over-plotted on a density map of our 
reference solution. The optical depth achieved in the wind is, at 
best, a few percent of that in the disc and therefore expect the 
disc ionisation to rely solely on its density distribution, regard- 
less of the properties of its wind. 

For a given accretion rate, the local disc density depends 
on the amplitude of the accretion speed Ur - nisCs - niseVk 
where m, is the sonic Mach number (measured here at the 
disc midplane), c, the sound speed and Va the Keplerian ve- 
locity. The main and major property of a JED is that accre- 
tion proceeds roughly at a sonic pace, namely m^ ~ 1 . This is 
a consequence of having magnetic fields close to equipartition 
(iFerreira & Pelletieilll995h and explains why cold and warm so- 
lutions give similar results. 

On the other hand, for a given accretion rate, the disc aspect 
ratio e is a critical parameter The density in the disc scales as e"^ 
(see Paper I) and the optical depth increases when e decreases, 
resulting in a lower disc midplane ionisation rate. We conclude 
that the disc thickness in any MAES solution is the key element 
to the disc ionisation and that warm or cold solutions with the 
same aspect ratio are readily interchangeable in that respect. 

4.4. JED ionisation during protostellar evolution 

So far, the paper has dealt with the prototype of a Class I YSO 
(M = 0.5 Mo, Mace = 10"^ Ma yr"') for which the JED cal- 
culation in paper I gave a disc aspect ratio e ~ 0.03. A warm 
MAES solution having both the right value of e and appropriate 
Class I/II jet properties has been used to derive the ionisation 
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Fig. 8. The disc aspect ratios obtained for parameters typical of 
Class 0, I and II protostars (see text for details), using the ap- 
proach developed in paper I. The curves end at the radius beyond 
which the JED becomes optically thin. 



degree in the JED of such a source. In order to extend this to 
the other stages of star formation, namely Class and Class II 
protostars we change M and Mace accordingly and calculate the 
new values of e following the method in paper I. As illustrated 
in fig.il 



- For a Class object with M - OA Mq and Mace 
1 0"^ Mq yr" ' , we find an average JED aspect ratio of e 

- For a Class II object with M = 0.5 Mq and Mace 
10"^ Mq yr"', we obtain e ~ 0.01. 



0.1. 



As discussed in 34.3l the disc ionisation state is independent of 
the jet properties. We therefore use isothermal (cold) MAES so- 
lutions with these values of e as they involve no extra physical 
ingredient (such as additional heat input at the disc surface). 

Figure |9] shows the disc midplane ionisation rate (left panel) 
and the ratio of the ionisation rate over its critical value (right 
panel) for Class O/I/II YSO prototypes, in solid, dashed and dot- 
ted lines respectively. The ionisation rate increases when moving 
from the youngest to the most evolved protostars (left). This is 
easily explained as Class discs are both denser and thicker than 
their more evolved counterparts, resulting in a smaller fraction of 
X-rays reaching the disc's midplane. 

The ratio (xl^^^ is shown on the right panel for both qui- 
escent (black) and flaring (blue) configuration of the X-ray 
source, calculated for = 1. The red dash-dotted lines cor- 
respond to the dead zone existence threshold for = 1 or 
- 100- Typical conservative quiescent X-ray properties 
(Lx = 10^^ erg s"', kTx - 3 keV, Zsi-c = 3 Rq) make the entire 
disc MRI-active if Kf^' = 1 but inactive if = 100, what- 
ever the stage of stellar formation; this is especially true for the 
Class objects. However and whatever the threshold, all types 
of protostars are MRI-active up to the disc's midplane when 
considering X-ray flaring characteristics (Lx - 10^'^ erg s"', 

^rx = 4keV,Zsi-c = 10 /^o). 

In the light of these results it appears that JEDs are very un- 
Ukely locations for dead zones to exist, whatever the evolution 
of the protostar may be. Class JEDs are more prone to host 
one than Class l/ll but this remains marginal and highly depen- 
dent on the quiescence X-ray source location and value of criti- 
cal magnetic number. 
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Fig. 9. Left: Log-ionisation rate (s ') on the disc midplane for a 
typical Class (solid line). Class I (dashed) and Class II (dotted) 
YSO. X-ray characteristic quiescence properties have been used 
{Lx - 10^^ erg s"', kTx - 3 keV, z^c - 3 Rq). Right: ionisation 
rate to its critical value for the same three star formation stages, 
calculated for - 1 ■ Both quiescence (black) and flaring X- 
rays (blue, Lx = 10^^" erg s"', kTx = 4 keV, z^i-c = 10 Rq) have 
been considered. The two dash-dotted red lines correspond to the 
threshold below which a dead zone is present, whether — 1 
orKT = 100. 



5. Discussion and concluding remarl^s 

We have revisited the question of the dead zones in accretion 
discs around young stellar objects by considering a new class 
of accretion solutions. Jet Emitting Discs. A JED is an accre- 
tion discs threaded by a large scale magnetic field, giving rise 
to self-confined jets carrying away a significant fraction of the 
disc angular momentum and energy. The jet emitting disc model 
relies on the self-similar MAES solutions, the magnetic fields of 
which are found smaller than, but close to, equipartition. This is 
close to the marginal stability criterion of the MRI, but still in 
the regime that allows it to be sustained. 

Only a small fraction of the released accretion energy of JED 
is converted into heat. For a given accretion rate, a JED is there- 
fore thinner than a standard accretion disc (SAD). Moreover, be- 
cause the torque due to the jets is dominant, accretion proceeds 
at a speed close to the speed of sound. This, in turn, leads to discs 
with a column density much smaller than in SADs. This directly 
translates into a more penetrating X-ray radiation, resulting in a 
higher ionisation degree compared to that of SADs. 

The ionisation structure of JEDs has been calculated using 
self-similar solutions along with the knowledge of their thermal 
structure. The radial extent of a JED is unknown as it depends 
on the importance of the magnetic flux available in the central 
regions of the disc. Thus, it certainly varies from one object to 
another. Moreover, for any given object, this extension will prob- 
ably also vary in time from Class to Class II. In this paper, the 
JED extent has been arbitrarily chosen between - 0.04 AU 
and o = 1 AU only, as it corresponds to the location where Class 
I/II jets seem to be launched: rj - 1-3 AU in the Class I pro- 
tostar DG Tau and rj ~ a few 0. 1 AU for some Class II YSO 
(lAnderson et al.ll2003l: iFerreira et al1l2006h . Wifli fliis caveat in 
mind and given the X-ray properties of young stellar objects, this 
work has mostly ruled out the existence of dead zones in JEDs. 

We did not include grain chemistry in our calculation argu- 
ing it may not be relevant to JEDs, as they occupy the innermost 
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regions of the discs. Previous studies have also shown gas-only 
and gas-grain chemical networks to converge in these innermost 
regions, comforting our assumption. However, the earliest stages 
of star formation may have JED extending up to a few AU, a 
distance at which grains could have a significant effect. Only a 
model using both gas-phase and gas-grain chemistry, be it a nu- 
merical chemical network calculation (e.g. Sem enov et al.l2004l) 
or a semi-analytical approach (iOkuzumi 2009) could settle this 
issue and this should be considered in future work. 

This result must also be modulated by our ignorance of the 
critical value of the magnetic Reynolds number above which the 
disc is MRl-active. Indeed, taking 7?™' = 1 allows the entire 
JED to be active, even during the - relatively - soft and low- 
luminosity quiescent state of YSOs. However, if "R™' = 100, 
then the critical ionisation degree is only achieved during X- 
ray flares. In that case, MRl can be continuously active in the 
entire JED only if the recombination timescale is longer than 
the quiescence phase between two flares, which typically lasts 
a few days. This is easily done for low values of the recombi- 
nation parameter JSq (such as in the metal dominated case where 
y6o = 3 X 10"" cm-' s"'), but becomes more problematic for high 
values (i.e., /3{) - 3 x 10"^ cm-' s"'). In this case, we cannot ex- 
clude the possibility for a dead zone to be intermittently present 
in JEDs during the X-ray protostellar quiescent phases and this 
effect should be stronger in the youngest (Class 0) objects. That 
matter requires a full time-dependent treatment that is beyond 
the scope of the present paper. 

The absence of a dead zone in the innermost regions of cir- 
cumstellar accretion discs depends on the outer extent ry of the 
JED. The direct observation of jet footpoints is not yet within 
the reach of current instruments. In the near future, ALMA 
(Atacama Large Millimeter Array) may be able to directly re- 
solve the launching regions of the closest and most powerful 
YSO jets, but in the meantime only indirect methods can be used. 
The value of ry and how it varies in time is an open issue as it 
depends on the magnetic field transport in discs. Since this trans- 
port is also dependent on the ionisation properties of the disc, 
we are facing here a typical "chicken or the egg" dilemma until 
a self-consistent treatment is made. For Class 0/1 objects, stan- 
dard approaches provide an inner dead zone boundary that can 
be significantly smaller than 1 AU, our typical ry. A SAD-JED 
radial transition could then be operated if the vertical magnetic 
field reaches a value close to equipartition. Within this frame- 
work, an outer dead zone would directly connect to an inner 
MRJ-active JED. This requires of course enough field advec- 
tion in the outermost SAD and through the dead zone. This pos- 
sibility clearly deserves fur ther investigations (see for instance 
iRothstein & Lovelacell2008h . 

The general trend is nevertheless that Jet Emitting Discs 
have smaller densities than SADs, sonic accretion and are com- 
pletely MRI-active: this makes them a hostile environment for 
planet formation. As a consequence, that process must be initi- 
ated and largely completed in the outer regions beyond the JED. 
Remarkably, several authors have pointed out the important role 
that the dead zone, and especially its innermost boundary, plays 
in both planet formation and migration processes. Indeed, solid 
materials may be trapped at its inner edge le ading to the for- 
mation of planet ary cores ( Kretke et al. 2009^, possibly within 
Rossby vortices (IVarniere & Tagger. ,2006) . This inner edge is 
systematically associated with a sudden decrease of the disc col- 
umn density and disc scale height. It thereby provides a very 
efficient means to stop the inward type I migration of proto- 
planetar y cores and may be a key ingredient to sav e planetary 
systems dMasset et al.ll2006t iMatsumura et al.ll2009l) . Note that 



these effects would be also at work at a SAD-JED transition. 
If planet formation starts early on in the protostellar evolution 
(Class O/I), when our results suggest that the dead zone should be 
truncated by a JED, then the conclusions of the above studies, if 
qualitatively unchanged, should be quantitatively reconsidered. 
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